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INTRODUCTION 

Natural  convection  and  radiation  modes  of  heat  transfer  are  commonly  applied  cooling 
techniques  for  electronic  equipment  of  low  to  moderate  power  density,  such  as  computer  chips, 
electronics,  and  telecommunication  boxes.  The  main  advantages  of  natural  convection  are  high 
reliability,  low  noise,  and  low  power  consumption.  Air  movement  in  a  naturally  ventilated  cabinet  is 
induced  by  a  pressure  difference  between  the  interior  and  exterior  of  the  cabinet.  This  pressure 
difference  occurs  due  to  a  density  gradient  in  the  air,  which  produces  a  force  that  pushes  the  air 
upward  and  out.  The  thermal  design  of  the  electronic  equipment  must  be  optimized  because  of  the 
relatively  low  efficiency.  To  generate  airflow  that  results  in  the  appropriate  temperature  rise  inside 
the  enclosure,  special  attention  must  be  paid  to  the  geometrical  configuration  of  the  heat  source  and 
the  ventilated  enclosure. 

Heat  sinks  are  devices  that  are  used  to  enhance  heat  dissipation  from  hot  surfaces  to  cooler 
air.  Typically,  the  fins  are  oriented  in  a  way  to  permit  a  natural  convection  air  draft  to  flow  upward 
through  the  rectangular  U-channels,  or  ducts,  formed  by  the  fins.  Heat  sink  design  goals  may  vary, 
but  in  this  particular  case,  optimization  of  the  heat  sink  is  the  main  concern.  One  of  the  most 
frequently  quoted  references  for  heat  sinks  is  Elenbass's  study  (ref.  1)  of  flow  between  parallel 
plates.  Elenbass  found  that  optimal  spacing  for  natural  convection  between  parallel  plates  was: 

Sop,  50 

H  Or  Pr  (1 

Considering  the  case  of  a  vertical  heat  sink  with  rectangular  fins,  equation  1  could  not  be 
used  without  modifications.  Additionally,  radiant  heat  transfer  also  has  to  be  considered  due  to  a  big 
part  of  heat  being  transferred  from  the  heat  sink  by  that  mode.  The  optimization  of  the  heat  sink  has 
become  an  essential  practice  for  optimum  performance.  In  present  cases,  optimization  of  the  heat 
sink  consisting  of  an  array  of  straight,  vertical,  and  rectangular  fins  is  performed  (fig.  1). 
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Figure  1 

Heat  sink  with  vertical  orientation 


MODEL  FORMULATION 

Heat  transfer  from  the  heat  sink  consists  of  radiation  and  convection  from  both  the  intra-fin 
passages  and  the  unshielded  surfaces  of  two  outer  fins.  In  this  report,  both  parts  are  considered 
separately. 

Intra-fin  Passages 

Convection 

For  the  U-channels  formed  by  rectangular  vertical  fins,  the  Nusselt  number  is 
calculated  with  a  modified  Elenbass's  equation  proposed  by  Van  De  Pol  and  Tierney  (ref.  2): 
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L 

y  =  -11.8  (in-1) 


All  physical  properties  except  p  are  evaluated  at  the  surface  temperature;  p  is  evaluated  at 
the  ambient  fluid  temperature. 

Once  the  Nusselt  number  is  obtained  using  equation  2,  the  convective  heat  transfer 
coefficient  for  the  U-channels  can  be  calculated  by: 


Nurk 

r 


Fin  efficiency  may  be  calculated  as: 


kHtf 


1 


2  hciHL- 


tanh 


2  hciHL~ 


kHtf 


Conductance  of  the  inter-fin  passages  now  can  be  calculated  as: 


cd  =  >AA 
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where  Ai  is  inter-fin  passage  area  given  as: 


A,  =  WH 


l  +  fcl)  L 
w 


(6) 


Radiation 

The  procedure  for  calculating  the  radiative  heat  transfer  coefficient  for  the  U-channels 
is  outlined  in  the  book  by  Ellison  (ref.  3).  It  is  assumed  that  all  surfaces  of  the  heat  sink  are  gray  and 
reflective  with  the  base  and  fins  at  a  uniform  temperature.  Figure  2  shows  a  single  U-channel  with 
numbered  interior  surfaces.  Numerals  1 , 3,  and  4  identify  the  heat  sink  surfaces  while  numerals  2,  5, 
and  6  refer  to  the  nonreflecting  ambient  air.  A  simplified  thermal  network  for  the  U-channels  is  shown 
in  figure  3. 


Note:  Numbers  are  centered  in  the  surface  that  they 
represent. 


Figure  2 

U-channel  interior  surface  identification 
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Figure  3 

Simplified  thermal  network  for  U-channel  radiation 
The  net  shape  factor,  F,  for  the  gray-body  U-channel  interior  is  given  as: 


F  = 


2  C 


NET 


[H(S  +  2L)] 


(7) 


where: 


C 


NET 


\iR.+Rt+RjR,+RJ+R,)-R,1 


(A  +M(A,  +  A  +*J(A  +  A,  +  *.)- A2  \- 

rMr,  +  Rd  +R,  )+RJtA-RjRjR.  +R„+R,  )+R„R,  ] 


is  the  overall  conductance  for  the  U-channel,  and  resistances  in  the  simplified  thermal  network  are 
given  in  the  following  form: 
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Rb 

K 

Rd 

Re 


1-6 

sA3 

2(1  ~g) 
sAl 

1 

AlFn  ^A3F35 
2 

A\F[i  +  2A[F15 

1 
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The  view  (shape  or  configuration)  factors  FI  2,  FI  3,  FI  5,  and  F35  are  computed  using  the 
appropriate  parallel  and  perpendicular  formulae  for  two  rectangles.  These  formulae  can  be  found  in 
any  heat  transfer  book  [e.g.,  A.F.  Mills  (ref.  4)]. 

Now,  using  the  result  obtained  by  equation  7,  the  radiation  conductance  can  be  calculated 
as: 


Cn  =  3.657*10“"  FA,. 


(r,+273)4  -(r„+273)4 

T,  -T. 


The  overall  conductance  for  the  inter-fin  passages  now  can  be  calculated  as: 


c  =c 

i  ci 


+  c 


ri 


(8) 


(9) 


Unshielded  Surfaces 

The  outer  surfaces  of  the  heat  sink  must  be  considered  in  order  to  obtain  more  realistic 
results.  Their  surface  conductance  is  typically  given  as: 


where: 


C  =  hA 


(10) 


and: 


Ko  + 


A0  =  2  HL 
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The  convective  heat  transfer  coefficient  for  the  vertical  fins  is  defined  as: 


rAT\0-25 


hco  =  0.0024 

j 

while  the  radiation  heat  transfer  coefficient  may  be  calculated  by: 

.„(r,+273)4-te+273)4 


h.  =  3.657x10 


T-T 


(11) 


(12) 


Overall  conductance  for  the  heat  sink  now  is: 

c,  =  Ct  +  c, 


(13) 


A  parametric  study  of  the  heat  sink  also  should  be  performed  to  determine  the  effects  of 
various  design  parameters  on  thermal  performance.  The  formulated  procedure  previously  shown  can 
be  coded  into  the  computer  program  using  various  programming  languages  such  as  Fortran,  C++, 
Pascal  or,  in  this  case,  it  can  be  done  with  Microsoft  Excel. 


MODEL  VALIDATION 

Validation  of  the  heat  sink  model  without  enclosure  was  done  using  the  following  parameters: 

•  Geometry 

•  Equipment  box 


• 

Height 

Hs  =  3.94  in. 

• 

Width 

Ws  =  5.42  in 

Fins 

• 

Height 

Hf  =  0.98  in. 

• 

Length 

Lf  =  3.94  in. 

Fin  thickness  and  conductivity  vary  depending  on  the  heat  sink  material  as  shown  in  table  1 . 

Table  1 

Varying  fin  thicknesses  and  conductivities  based  on  material 


Fin  material 

Thickness  of  the  fin,  tf  (in.) 

Fin  conductivity,  k(W/in.K) 

Aluminum 

0.021 

4.84 

Steel 

0.046 

1.625 

Porcelain 

0.094 

0.267 

•  Ambient  temperature:  Ta=  20°C 

•  Fin  temperature:  Tf=50°C 

•  Fin  emissivity:  s  =0.82 

Table  2  gives  a  comparison  of  the  number  of  fins  and  heat  dissipation  by  fin  between  the 
present  analysis  and  the  results  obtained  from  Aihara  (ref.  5). 
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Table  2 

Comparison  of  results  for  optimal  number  of  fins 


Materia  of  No.  of  fins  Q  dissipation/fin, W 


heat  sink 

Aiharas 

Present 

Aiharas 

Present 

Aluminum 

19 

18 

17.9 

17.7 

Steel 

17 

16 

16.8 

16.6 

Porcelain 

15 

14 

15.1 

14.6 

The  results  of  table  2  agree  with  the  results  from  Aiharas,  and  the  differences  are  actually 
due  to  the  efficiencies  being  accounted  for  in  this  present  calculation. 


PROBLEM  STATEMENT 

An  electronic  box  of  the  dimensions  24  in.  by  18  in.  by  10  in.  (H  x  W  x  D)  dissipates  250  W  of 
heat.  Natural  convection  and  radiation  are  used  to  cool  the  external  surfaces  of  the  box.  The 
maximum  expected  ambient  temperature  is  46°C,  and  the  maximum  allowable  box  surface 
temperature  is  81  °C.  It  is  assumed  that  all  heat  is  dissipated  only  through  the  wall  with  dimensions  of 
24  in.  by  18  in.  A  heat  sink  with  rectangular  fins  of  dimensions  24  in.  by  1.77  in.  by  0.12  in.  [H  x  L  x 
fin  thickness  (tf)]  is  attached  to  the  box.  The  fin  material  is  aluminum  (k=4.84  W/in.K)  with  fin 
emissivity  s=0.9.  The  objective  is  to  optimize  the  rectangular  heat  sink  to  meet  design  requirements. 
Additionally,  a  parameterization  study  should  be  performed  to  assess  heat  sink  efficiency. 


RESULTS 

Figure  4  depicts  the  average  heat  sink  heat  dissipation  and  conductance  as  a  function  of 
spacing.  While  the  fin  spacing  is  large,  the  conductance  of  the  heat  sink  stays  relatively  constant.  By 
decreasing  the  spacing  between  the  fins,  more  fins  can  be  added  to  the  base  plate,  thereby 
increasing  the  total  surface  area  for  heat  transfer.  Since  the  total  heat  flow  is  proportional  to 
conductance  and  total  surface  area,  decreasing  the  spacing  further  still  improves  total  heat  transfer 
until  optimum  spacing  is  reached.  At  the  point  of  optimum  spacing,  heat  transfer  from  the  sink 
reaches  its  maximum  value.  Below  optimum  spacing,  the  total  heat  transfer  falls  rapidly  due  to  fin 
interference  and  choked  airflow. 
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Figure  4 

Heat  sink  heat  dissipation  and  conductance  as  a  function  of  spacing 

The  conductance  of  heat  sink  intra-fin  passages  is  shown  in  figure  5  for  different  numbers  of 
fins.  By  increasing  the  number  of  fins,  conductance  due  to  convection  (cci)  steadily  increases  until 
an  optimum  number  of  fins  is  reached.  Beyond  that  number,  the  convection  conductance  falls 
steadily.  Conductance  due  to  radiation  (cri)  declines  with  an  increase  in  the  number  of  fins  as  well. 
Overall,  intra-fin  conductance  is  obtained  by  adding  conductance  due  to  both  radiation  and 
convection. 


ci 

cri 


cci 


Number  of  fins 


Figure  5 

Conductance  of  intra-fin  passages  as  a  function  of  the  fin  number 
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Figure  6  shows  the  dependence  of  the  heat  sink's  intra-fin  passages  conductance  on 
spacing.  As  shown,  large  spacing  conductance  stays  relatively  constant  until  a  certain  limit  is 
reached.  By  continuing  to  decrease  the  fin  spacing,  conductance  starts  to  increase  rapidly  until  the 
optimum  spacing  is  reached.  At  this  point,  the  highest  conductance  and  heat  dissipation  from  the 
heat  sink  is  achieved.  Decreasing  the  fin  spacing  further  causes  a  sudden  fall  in  both  conductance 
and  heat  dissipation,  and  the  advantages  of  the  heat  sink  loses  its  significance. 


Considering  total  heat  sink  conductance  and  heat  transfer  (fig.  7),  the  optimal  number  of  fins 
is  28  with  a  maximum  heat  transfer  of  298.28  W  and  conductance  of  8.52  W/K. 


Number  of  fins  N 


Figure  7 

Heat  transfer  and  conductance  of  the  heat  sink 
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Figure  8  shows  the  lowest  temperature  of  ~75°C  is  achieved  with  28  fins,  which  is  well  below 
the  design  requirement  of  81  °C. 


Number  of  the  fins 


Figure  8 

Temperature  at  the  base  of  the  heat  sink 

Figure  9  represents  the  dependence  of  the  heat  sink  base  temperature  on  spacing.  Again,  as 
the  fin  spacing  is  reduced,  the  base  temperature  changes  slowly  at  first  and  ramps  rapidly.  This 
trend  reaches  its  limit  at  the  point  of  optimal  spacing  when  the  minimum  temperature  of  the  heat  sink 
base  occurs.  If  the  decrease  in  spacing  is  continued,  sudden  increase  in  temperature  of  the  base 
occurs  due  to  choked  airflow  and  fin  interference.  The  lowest  temperature  appears  at  the  spacing  of 
0.5422  in. 
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spacing  (in) 


Figure  9 

Temperature  at  the  base  of  the  heat  sink  as  a  function  of  fin  spacing 

Parametric  Study 

An  extensive  parametric  study  of  the  heat  sink  is  now  performed  to  investigate  possible 
material  and  space  savings  in  the  design  process.  The  first  parameter  considered  was  the  thickness 
of  the  fin.  Starting  with  a  small  number  of  fins,  overall  conductance  of  the  heat  sink  was  the  same 
regardless  of  fin  thickness.  By  increasing  the  number  of  fins,  the  conductance  of  the  heat  sink 
becomes  larger  for  thinner  fins  (fig.  1 0).  Also,  the  optimal  number  of  fins  increases  when  compared 
to  the  thicker  fins.  The  temperature  of  the  base  decreased  with  a  decrease  in  fin  thickness  (fig.  11), 
and  the  heat  transfer  from  the  heat  sink  increased  with  the  decreased  thickness  (fig.  12).  From  these 
findings,  it  can  be  concluded  that  thinner  fins  transferred  more  heat  due  to  the  increased  number  of 
fins. 
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Figure  10 

Conductance  of  the  heat  sink  for  different  fin  thicknesses  (tf) 
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Figure  1 1 

Temperature  at  the  base  of  the  heat  sink  for  different  fin  thickness  (tf) 
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Figure  12 

Heat  transfer  from  the  heat  sink  for  different  fin  thicknesses  (tf) 

The  next  parameter  considered  was  fin  height,  Hf.  Smaller  height  has  a  negative  influence 
overall  heat  sink  conductance  (fig.  13)  causing  an  increase  in  fin  base  temperature  (fig.  14)  and 
decreases  the  overall  heat  transfer  from  the  sink  (fig.  15).  It  is  impossible  to  meet  the  design 
temperature  at  the  base  of  the  heat  sink  by  decreasing  the  height  of  the  fin  too  much  (in  this  case, 
decreasing  Hf  below  18  in.). 
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Figure  13 

Conductance  of  the  heat  sink  for  different  fin  height  (L) 
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Figure  14 

Temperature  at  the  base  of  the  heat  sink  for  different  fin  height  (L) 
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Figure  15 

Heat  transfer  from  the  heat  sink  for  different  fin  height  (L) 

The  last  parameter  to  consider  was  fin  length,  L.  Fin  length  has  an  even  more  negative 
influence  on  the  behavior  of  the  heat  sink  than  the  height.  Decreasing  the  length  of  a  fin  significantly 
reduces  overall  heat  sink  conductance  (fig.  16),  and  as  a  result,  the  temperature  of  the  base 
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increases  above  the  required  design  temperature  of  81  °C  (fig.  17).  The  heat  sink’s  ability  to  dissipate 
heat  from  the  chip  or  electronics  is  also  greatly  reduced  (fig.  18). 
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Figure  16 

Conductance  of  the  heat  sink  for  different  fin  length  (H) 


Number  of  fins 


L=1 .77" 
L=1 .38" 
L=1 " 


Figure  17 

Temperature  at  the  base  of  the  heat  sink  for  different  fin  length  (H) 
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Figure  18 

Heat  transfer  from  the  heat  sink  for  different  fin  length  (H) 


CONCLUSIONS 

Heat  sink  optimization  helps  us  to  determine  the  dimensions,  weight,  and  thermal 
performance  in  order  to  meet  design  requirements.  As  shown,  finding  the  optimal  number  of  fins 
determines  the  best  solution  for  temperature  at  the  base  of  the  sink.  Also,  it  became  clear  that  when 
the  number  of  fins  was  increased  past  an  optimal  value,  the  thermal  performance  of  the  sink  would 
worsen,  leading  to  choked  flow  inside  the  U-channels  and,  eventually,  electronic  box  failure.  The 
goal  is  to  determine  the  optimal  number  of  fins  in  order  to  maximize  heat  transfer.  Thickness  of  the 
fin  is  the  major  parameter  to  heat  sink  performance.  Although  thinner  fins  increase  the  optimum 
number  of  the  fins,  they  actually  decrease  the  weight  of  the  heat  sink  and  improve  the  heat  sink’s 
thermal  performance. 

It  can  be  concluded  that  although  all  assumptions  (uniform  temperature  at  the  base  and  the 
fins,  gray  surfaces,  etc.)  made  in  the  model  formulation  are  not  satisfied  perfectly,  they  are  more 
than  satisfactory  for  most  engineering  calculations. 
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GLOSSARY 

(B  Thermal  coefficient  of  volume  expansion 

c  Emissivity 

F  Shape  factor 

Gr  Grashof  number 

H  Fin  height,  (in.) 

k  Thermal  conductivity,  (W/in.-K) 

r|  Fin  efficiency 

Nu  Nusselt  number 

Pr  Prandtl  number 

Q  Heat  dissipation,  (W) 

Ra  Rayleigh  number 

R  Thermal  resistances,  (K/W) 

S  Fin  spacing,  (in.) 

tf  Fin  thickness,  (in.) 

Ta  Ambient  temperature,  (K) 

Ts  Surface  temperature,  (K) 

Hf  Fin  height 
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